Absolute Free Energy of Binding Calculations for Macrophage Migration Inhibitory Factor in Complex with a Druglike Inhibitor

Calculation of the absolute free energy of binding (ΔGbind) for a complex in solution is challenging owing to the need for adequate configurational sampling and an accurate energetic description, typically with a force field (FF). In this study, Monte Carlo (MC) simulations with improved side-chain and backbone sampling are used to assess ΔGbind for the complex of a druglike inhibitor (MIF180) with the protein macrophage migration inhibitory factor (MIF) using free energy perturbation (FEP) calculations. For comparison, molecular dynamics (MD) simulations were employed as an alternative sampling method for the same system. With the OPLS-AA/M FF and CM5 atomic charges for the inhibitor, the ΔGbind results from the MC/FEP and MD/FEP simulations, −8.80 ± 0.74 and −8.46 ± 0.85 kcal/mol, agree well with each other and with the experimental value of −8.98 ± 0.28 kcal/mol. The convergence of the results and analysis of the trajectories indicate that sufficient sampling was achieved for both approaches. Repeating the MD/FEP calculations using current versions of the CHARMM and AMBER FFs led to a 6 kcal/mol range of computed ΔGbind. These results show that calculation of accurate ΔGbind for large ligands is both feasible and numerically equivalent, within error limits, using either methodology

Complementarity Between a Docking and a High-Throughput Screen in Discovering New Cruzain Inhibitors†

Virtual and high-throughput screens (HTS) should have complementary strengths and weaknesses, but studies that prospectively and comprehensively compare them are rare. We undertook a parallel docking and HTS screen of 197861 compounds against cruzain, a thiol protease target for Chagas disease, looking for reversible, competitive inhibitors. On workup, 99 % of the hits were eliminated as false positives, yielding 146 well-behaved, competitive ligands. These fell into five chemotypes: two were prioritized by scoring among the top 0.1 % of the docking-ranked library, two were prioritized by behavior in the HTS and by clustering, and one chemotype was prioritized by both approaches. Determination of an inhibitor/cruzain crystal structure and comparison of the high-scoring docking hits to experiment illuminated the origins of docking false-negatives and false-positives. Prioritizing molecules that are both predicted by docking and are HTS-active yields well-behaved molecules, relatively unobscured by the false-positives to which both techniques are individually prone

Stability mechanisms of a thermophilic laccase probed by molecular dynamics.

Laccases are highly stable, industrially important enzymes capable of oxidizing a large range of substrates. Causes for their stability are, as for other proteins, poorly understood. In this work, multiple-seed molecular dynamics (MD) was applied to a Trametes versicolor laccase in response to variable ionic strengths, temperatures, and glycosylation status. Near-physiological conditions provided excellent agreement with the crystal structure (average RMSD ∼0.92 Å) and residual agreement with experimental B-factors. The persistence of backbone hydrogen bonds was identified as a key descriptor of structural response to environment, whereas solvent-accessibility, radius of gyration, and fluctuations were only locally relevant. Backbone hydrogen bonds decreased systematically with temperature in all simulations (∼9 per 50 K), probing structural changes associated with enthalpy-entropy compensation. Approaching T opt (∼350 K) from 300 K, this change correlated with a beginning "unzipping" of critical β-sheets. 0 M ionic strength triggered partial denucleation of the C-terminal (known experimentally to be sensitive) at 400 K, suggesting a general salt stabilization effect. In contrast, F(-) (but not Cl(-)) specifically impaired secondary structure by formation of strong hydrogen bonds with backbone NH, providing a mechanism for experimentally observed small anion destabilization, potentially remedied by site-directed mutagenesis at critical intrusion sites. N-glycosylation was found to support structural integrity by increasing persistent backbone hydrogen bonds by ∼4 across simulations, mainly via prevention of F(-) intrusion. Hydrogen-bond loss in distinct loop regions and ends of critical β-sheets suggest potential strategies for laboratory optimization of these industrially important enzymes

In Silico Improvement of beta(3)-Peptide Inhibitors of p53 center dot hDM2 and p53 center dot hDMX

There is great interest in molecules capable of inhibiting the interactions between p53 and its negative regulators hDM2 and hDMX, as these molecules have validated potential against cancers in which one or both oncoproteins are overexpressed. We reported previously that appropriately substituted β(3)-peptides inhibit these interactions and, more recently, that minimally cationic β(3)-peptides are sufficiently cell permeable to upregulate p53-dependent genes in live cells. These observations, coupled with the known stability of β-peptides in a cellular environment, and the recently reported structures of hDM2 and hDMX, motivated us to exploit computational modeling to identify β-peptides with improved potency and/or selectivity. This exercise successfully identified a new β(3)-peptide, β53-16, that possesses the highly desirable attribute of high affinity for both hDM2 as well as hDMX and identifies the 3,4-dichlorophenyl moiety as a novel determinant of hDMX affinity. [Image: see text

Unfolding of the Amyloid β-Peptide Central Helix: Mechanistic Insights from Molecular Dynamics Simulations

Alzheimer's disease (AD) pathogenesis is associated with formation of amyloid fibrils caused by polymerization of the amyloid β-peptide (Aβ), which is a process that requires unfolding of the native helical structure of Aβ. According to recent experimental studies, stabilization of the Aβ central helix is effective in preventing Aβ polymerization into toxic assemblies. To uncover the fundamental mechanism of unfolding of the Aβ central helix, we performed molecular dynamics simulations for wild-type (WT), V18A/F19A/F20A mutant (MA), and V18L/F19L/F20L mutant (ML) models of the Aβ central helix. It was quantitatively demonstrated that the stability of the α-helical conformation of both MA and ML is higher than that of WT, indicating that the α-helical propensity of the three nonpolar residues (18, 19, and 20) is the main factor for the stability of the whole Aβ central helix and that their hydrophobicity plays a secondary role. WT was found to completely unfold by a three-step mechanism: 1) loss of α-helical backbone hydrogen bonds, 2) strong interactions between nonpolar sidechains, and 3) strong interactions between polar sidechains. WT did not completely unfold in cases when any of the three steps was omitted. MA and ML did not completely unfold mainly due to the lack of the first step. This suggests that disturbances in any of the three steps would be effective in inhibiting the unfolding of the Aβ central helix. Our findings would pave the way for design of new drugs to prevent or retard AD

Quantifying Water-Mediated Protein–Ligand Interactions in a Glutamate Receptor: A DFT Study

It is becoming increasingly clear that careful treatment of water molecules in ligand–protein interactions is required in many cases if the correct binding pose is to be identified in molecular docking. Water can form complex bridging networks and can play a critical role in dictating the binding mode of ligands. A particularly striking example of this can be found in the ionotropic glutamate receptors. Despite possessing similar chemical moieties, crystal structures of glutamate and α-amino-3-hydroxy-5-methyl-4-isoxazole-propionic acid (AMPA) in complex with the ligand-binding core of the GluA2 ionotropic glutamate receptor revealed, contrary to all expectation, two distinct modes of binding. The difference appears to be related to the position of water molecules within the binding pocket. However, it is unclear exactly what governs the preference for water molecules to occupy a particular site in any one binding mode. In this work we use density functional theory (DFT) calculations to investigate the interaction energies and polarization effects of the various components of the binding pocket. Our results show (i) the energetics of a key water molecule are more favorable for the site found in the glutamate-bound mode compared to the alternative site observed in the AMPA-bound mode, (ii) polarization effects are important for glutamate but less so for AMPA, (iii) ligand–system interaction energies alone can predict the correct binding mode for glutamate, but for AMPA alternative modes of binding have similar interaction energies, and (iv) the internal energy is a significant factor for AMPA but not for glutamate. We discuss the results within the broader context of rational drug-design

Conformational Changes and Slow Dynamics through Microsecond Polarized Atomistic Molecular Simulation of an Integral Kv1.2 Ion Channel

Structure and dynamics of voltage-gated ion channels, in particular the motion of the S4 helix, is a highly interesting and hotly debated topic in current membrane protein research. It has critical implications for insertion and stabilization of membrane proteins as well as for finding how transitions occur in membrane proteins—not to mention numerous applications in drug design. Here, we present a full 1 µs atomic-detail molecular dynamics simulation of an integral Kv1.2 ion channel, comprising 120,000 atoms. By applying 0.052 V/nm of hyperpolarization, we observe structural rearrangements, including up to 120° rotation of the S4 segment, changes in hydrogen-bonding patterns, but only low amounts of translation. A smaller rotation (∼35°) of the extracellular end of all S4 segments is present also in a reference 0.5 µs simulation without applied field, which indicates that the crystal structure might be slightly different from the natural state of the voltage sensor. The conformation change upon hyperpolarization is closely coupled to an increase in 310 helix contents in S4, starting from the intracellular side. This could support a model for transition from the crystal structure where the hyperpolarization destabilizes S4–lipid hydrogen bonds, which leads to the helix rotating to keep the arginine side chains away from the hydrophobic phase, and the driving force for final relaxation by downward translation is partly entropic, which would explain the slow process. The coordinates of the transmembrane part of the simulated channel actually stay closer to the recently determined higher-resolution Kv1.2 chimera channel than the starting structure for the entire second half of the simulation (0.5–1 µs). Together with lipids binding in matching positions and significant thinning of the membrane also observed in experiments, this provides additional support for the predictive power of microsecond-scale membrane protein simulations

Unfolding Simulations of Holomyoglobin from Four Mammals: Identification of Intermediates and β-Sheet Formation from Partially Unfolded States

Myoglobin (Mb) is a centrally important, widely studied mammalian protein. While much work has investigated multi-step unfolding of apoMb using acid or denaturant, holomyoglobin unfolding is poorly understood despite its biological relevance. We present here the first systematic unfolding simulations of holoMb and the first comparative study of unfolding of protein orthologs from different species (sperm whale, pig, horse, and harbor seal). We also provide new interpretations of experimental mean molecular ellipticities of myoglobin intermediates, notably correcting for random coil and number of helices in intermediates. The simulated holoproteins at 310 K displayed structures and dynamics in agreement with crystal structures (R g ~1.48-1.51 nm, helicity ~75%). At 400 K, heme was not lost, but some helix loss was observed in pig and horse, suggesting that these helices are less stable in terrestrial species. At 500 K, heme was lost within 1.0-3.7 ns. All four proteins displayed exponentially decaying helix structure within 20 ns. The C- and F-helices were lost quickly in all cases. Heme delayed helix loss, and sperm whale myoglobin exhibited highest retention of heme and D/E helices. Persistence of conformation (RMSD), secondary structure, and ellipticity between 2-11 ns was interpreted as intermediates of holoMb unfolding in all four species. The intermediates resemble those of apoMb notably in A and H helices, but differ substantially in the D-, E- and F-helices, which interact with heme. The identified mechanisms cast light on the role of metal/cofactor in poorly understood holoMb unfolding. We also observed β-sheet formation of several myoglobins at 500 K as seen experimentally, occurring after disruption of helices to a partially unfolded, globally disordered state; heme reduced this tendency and sperm-whale did not display any sheet propensity during the simulations